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Abstract 

A simple modification of the standard Renormalization Group (RG) technique for 
the study of quantum spin systems is introduced. Our method which takes into ac- 
count the effect of boundary conditions by employing the concept of superblock, may 
be regarded as a simple way for obtaining first estimates of many properties of spin 
systems. By applying this method to the xxz spin | Heisenberg chain, we obtain the 
ground state energy with much higher accuracy than the standard RG. We have also 
obtained the staggered magnetization and the z-component of spin-spin correlation 
function which confirms the absence of long-range order in the massless region of the 
ID xxz model. 
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1 Introduction 



Soon after the introduction of real space Renormalization Group (RG) by Wilson and it's 
application to the Kondo problem 1 , it was found that this type of RG does not always 
produce accurate numerical results, as compared with other powerful techniques. One of the 
weak performances of the RG method was found in Lee's work on 2D Anderson localization 2 . 
Using numerical RG method he concluded that there is a critical amplitude for the random 
potential which induces localization. However, not long after Lee's work, Lee and Fisher 2 , 
using a different approach, found that the 2D model is logarithmically localized even for 
arbitrary small randomness. This result is now generally accepted. Such unsatisfactory 
behavior of the RG method compared with the other methods such as Quantum Monte 
Carlo, was one of the reasons that the RG method for quantum lattice problems remained 
undeveloped during the 80 's. 

Let us first describe the standard RG method 3 ( see ref.4 for a modern treatment and 
references therein). First, one divides the whole lattice into small blocks (fig.l) and obtains 
the lowest states of each isolated block for a particular Boundary Condition (BC).The 
effect of inter-block interactions is taken into account by constructing an effective Hamilto- 
nian H e ^ which now acts on a smaller Hilbert space 7i e ^ embedded in the original one. In 
this new Hilbert space, each of the former blocks is treated as a single site . 
The technical way of implementing this idea is to construct an embedding operator 

T : H eff — ► H 

and a truncation operator 

T f : H — >H eff 
and demand the commutativity of the following diagram 4 : 

H eff n 
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i.e. : TH e ff = HT. From the last relation one obtains the effective Hamiltonian as : 

H eff = T ^ HT 

Note that the operators T and satisfy the relation T^T = l^e// but TT^ ^ 1^. More 
precisely one divides H as H = H B + H BB (fig.l) where H B = J2i hf is the sum of block 
Hamiltonians and H BB = J2i,j hf j is the sum of inter-block Hamiltonians and constructs 
T — YliTi where Tj is constructed as follows : 

m 

T> = ElO<i| (2) 
i=i 

where m is the number of low energy states kept. Note that each Tj acts trivially on all the 
other blocks J ^ I consequently [T I: Tj] = 0. Repeating the RG steps one hopes to restrict 
himself to spaces of lower and lower energy and finally arrives at the ground state energy. 
In the general case this process induces a renormalization of coupling constants : 

H e " = H(k' v k' 2 , ...) = T^H(k 1 , k 2 , ...)T (3) 

where K = (ki, k 2 , ■■■) is the coupling constant space. 

In principle the form the Hamiltonian changes in a such RG program rendering analytical 
treatment, rather difficult if not possible. In the standard RG treatment of some models 
including xxz Heisenberg chain , it is possible to choose the embedding operator so that no 
new terms appear in the renormalized Hamiltonian 3 . In such cases, repeating the RG steps, 
leads to the following general form : 

n-l aj 

T^" H(K)T n = £ — ef + H„(K™) (4) 

here H N is the Hamiltonian of the N-sites lattice, / is the number of sites in each block, n is 
the number of RG steps performed, e B (K^) is the ground state energy of each block after 
the p-th step and under RG : — > i^( n + 1 ) ) then the ground state of the whole lattice is 
then given by : 

OO AT 

£ = Ew^)) (5) 

p=0 1 

However the main difficulty of the method is that by fixing a particular BC on an isolated 
block one may loose a number of states which contribute to the ground state of the whole 
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lattice due to the interaction of the block with its surrounding. Stated in another way, the 
ground state wave function of the whole lattice for a particular BC is not a simple juxtapo- 
sition of the ground state wave functions of the blocks. This point is clearly highlighted in 
ID Tight Binding model in which the standard RG fails for some type of BCs 5 ' 6 . It must be 
noticed that this difficulty is not removed by increasing the size of blocks. 

In the Density Matrix Renormalization Group 7 one embeds each block into a larger(super) 
block and considers the block as a quantum system in interaction with a reservoir(the rest 
of the superblock). The block can then be described by a reduced density matrix, whose 
eigenkets with the highest eigenvalues are used to construct the embedding operator(T). 
Provided that one keeps a large number of states in each step of truncation(i.e. m = 16, 24, 36 
and 44 in four different runs 7 ), this method leads to very accurate results for energies and 
correlation functions of the isotropic xxx Heisenberg chain. 

However keeping a large number of states renders the problem of identification of coupling 
constants of the effective Hamiltonian very difficult, which means that one should make a 
compromise in obtaining accurate results for energies and correlation functions on the one 
hand and coupling constant trasformations on the other, the later objective is of course 
important if one is interested in studying the critical properties of these systems. In fact while 
the DMRG method is quite appropriate for obtaining the former properties, the standard 
RG appeares to give better results for the later properties 8 . One should adapt one or another 
of these RG schemes to obtain reasonable results in both directions. In ref.[9] a strategy has 
been suggested for the adaptation of DMRG so that one can also derive the transformation 
of coupling constants and the critical properties of the xxz Heisenberg chain. 

The present paper has a different starting point in that it adapts the standard RG and 
incorporates in it the interaction of the block with its environment, to do the same general 
task. The method presented in this paper and in ref. [9] are to be regarded as complementary. 
Depending on the problem and the desired accuracy either of them can be used. The 
simplicity of our method is such that one may hope to apply it to two dimensional systems. 
Although it does not yield as much accuracy as in DMRG, is much simpler practically, so that 
for obtaining first estimates of many properties of spin or fermion systems one may try this 
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method 10 . For the xxz chain, with blocks of 3 sites in the infinite chain and with retaining 
only m = 2 states in each block, we have obtained the ground state energy with an accuracy 
much better than the standard RG method. The results are compared in table-1 and table-2. 
We have also obtained the staggered magnetization and spin-spin correlation functions in 
the massless region of the xxz model. These results verify the absence of long-range order 
in this region which is in agreement with the known results. 

2 A modified RG for the s = \ xxz chain 

The system under consideration is the anti-ferromagnetic s — \ xxz chain with the following 
Hamiltonian : 

N-l 

H N (J, A) = (J/4) £ KX +1 + afaf +1 + A(a*a* +1 )) (6) 
i=i 

where J is positive and A is the anisotropy parameter which is limited by (0 < A < 1) in 
the massless region. 

Fig. 2 shows the decomposition of the chain into isolated blocks and super blocks with 
an odd number of sites per each block. Thus the block Hamiltonian h B and the superblock 
Hamiltonian h SB are 

h B = (J/4)(<j 2 X + of at + + + A(<7 2 V* + aloft) (7) 

h SB = (J/4) {alal + a\a\ + Aafa z 2 ) + h B + (J/4) {a%a% + a\a\ + Aa 4 V 5 2 ) (8) 

The basic steps of our method are follows : 

step-1) First one diagonalizes the superblock Hamiltonian (h SB ). Although h SB does 
not have su(2) symmetry, it has a number of symmetries which helps in it's diagonalization. 
Defining the z-component of total spin S z = \{a\ + °"2 + °3 + a l + a l)i the flip operator 
a = afa^a^afa^, and the parity operator(n), Il\aia2cx^a 4 a^) = \a 5 a4a 3 a 2 ai) . It is easy to 
check that, H commutes with n, a and S z . Furthermore [n, S z ] = 0, [n, a] = 0, which means 
that the energy eigenstates of h SB can be chosen to have definite parity and the z-component 
of spin. The ground state is doubly degenerate, one with the z-component of the spin equal 
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to +| and the other equal to — \. 

h SB \±b =efV,A)|±i> (9) 

2 SB 2 SB 

where | ± \) SB are the doubly degenerate ground state of h SB and ef B is their corresponding 
energy. We note that the states of lowest energy are even parity states with their z-component 



of spin equal to ±|. Thus the generic form of | + \) qR is as follows 

i) 

Tsb 



_ 2 . j-j-iiAo 6"' vl ' l/ I 1 2 /SB 

+ -) = a x | + + + ) +a 2 | + + --+) + a 3 | + + -+-) 



+ a 4 | + - + +-) + a 5 | + - + -+) + a 6 \ - + + +-) + a 4 | - + + -+) 

+ a 3 \ - + -++) + a 2 | + --++) + ai| -- + ++) (10) 

The state | — \) SB is obtained from | + |) by flipping all the spins. All the coefficient a, ; 
i — 1, ... ,6 can be expressed as functions of J and A although we do not need to explicitly 
write them at this step, and ef B , itself can be calculated numerically to any desired accuracy. 

step-2) Knowing that the ground state of the block is a 2-dimensional subspace with 
s z — ±| , we then project the | ± \) SB states onto the s z = ±| subspaces of the block 
Hilbert space. It is this step which effectively smoothes out the sharp effects of the boundary 
conditions, by immersing the block into a superblock (which to some extent simulates the 
effect of the rest of the lattice), the projected state does not restrict to any particular 
boundary conditions of an isolated block (open, periodic, ...). Moreover, in this way the 
form of the Hamiltonian does not change, as we will see later. What we do is that, after 
projection onto the block Hilbert space we implement the projection operator P + which is 
shown below, to reach the s z = +| subspaces of the block Hilbert space. 

P + = I + +-)(- + +1 + I + -+)(+- +| + | - ++)(+ + -| (11) 
Then the resulting normalized states are 

\ + h =P + \ + h =(l/A)((a 1 + a A )\ + +-)+2a 3 \ + -+) + (a 1 + a 4 )\-++)) 

2 B 2 SB 

\-h =P'\~h =(l/A)(( ai + a 4 )\- ■-+) + 2a 3 |-+-) + (ai + a 4 )| + — )) (12) 
2 b 2 sb 

where A is a normalization factor. Straightforward calculations give the following results 

o-n _ (s-A)(2e 2 +sA-2) / -i q\ 

P = 2a 3 - e(e 2_ A 2_ 1) — 
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where 

e = (2/J)ef . (14) 
Numerical calculation gives the following result for ef B (0 < A < 1), 

e s Q B (J, A) = J{a + 6A + cA 2 + dA 3 + eA 4 + /A 5 + 5A 6 + /iA 7 + kA 8 + /A 9 ) (15) 

with 

a = -1.36623 b = -0.63707 c = 2.92264 d = -22.07630 e = 84.52129 
/ = -187.92379 g = 251.64014 h = -200.06805 k = 86.98911 / = -15.92981. 

step-3) The embedding operator is now constructed in the following form 

Ti = (\ + h (+\) + (\-h (-D as) 

A B I B 

where |+) and |— ) are the renamed base kets of the effective Hamiltonian Hilbert space of 
each block. 

Having the form of the embedding operator, one can calculate the projection of any 
operator onto the effective Hilbert space, as in eq.(l) i.e. 

^eff = TVfT = ( 2^ )a V i = lt3;lM = Xiy 

or ,f = Tta*T = (&)a? < = 1,3 (17) 

where a\ acts as pauli matrices in the new Hilbert space. Fig.l shows that the interaction 
between blocks are in the following form 

hf B j = W)(ofrfj + ahah + Aa^alj). (18) 

Then the effective Hamiltonian between the new sites are 

h e /J = TjTjhf^TrTj (19) 

where I and J are two neighbouring blocks while Tj and Tj are their corresponding embedding 
operators. Doing so we arrive at : 

h e /J = ( J'/4) {of a] + a'?*'} + A'a'fa'f) (20) 



where 

A' = (^) 2 A (21) 

Eq.(21) describes the RG flow of J and A. The RG flow goes to the origin of J — A plane 
which verifies that we are in the massless region. 

step-4) The crudest estimate for the ground state energy is just the sum of ground state 
energies of all the blocks or superblocks. Thus 

(|)(°) = lef(J,A) (22) 

or 

(f) (0) = ^f(J,A) (23) 

where (^) < -°- ) is the ground state energy of the whole system in the zeroth order approxima- 
tion. Eq.(23) in fact defines ef as : 

ef = \4 B (24) 
However this estimate clearly neglects the contribution to the energy from interaction 
between blocks. To take this missing contribution into account, we add to the ground 

state energy of the new lattice, and obtain : 

(f)« = i e ?(J,A) + i e ?(J«,A( 1 )) (25) 
iterating this procedure we finally obtain 

JV n=0 ° 

or 

§ = t ^T(|«fV"',AW)) (27) 

where = J , A^ ^ = A and J*™) , A^ n ^ are the renormalized coupling constants after 
n-steps of RG. In fact the above argument describes the physical interpretation of the formal 
limiting process defined in ref.4, in which the ground state energy is obtained as : 

E = lim n ^T^ n HT n 



However in standard RG the terms e B (J : A) is the result of diagonalization of h B , i.e: 
T^h B T = e B , while in our method e B is just another expression for ^Cq 3 where e^ 3 is the 
eigenvalue of h SB , I 1 and I are the sizes of the superblocks and blocks respectively. The results 
for the ground state energy per site, obtained by this modified RG method are collected in 
table- 1, for a series of anisotropy parameter (0 < A < 1). The results are compared both 
with the exact results (obtained by the Bethe ansatz for xxx 11 and xxz 12 spin \ chains) and 
with those obtained from the standard RG 13 . It is clearly seen that this type of modified 
RG method yields much better results than the standard RG method. 

A remark is in order here, concerning the non-variational nature of our results, which is 
caused by our crude estimate for the ground state energy of each block in terms of that the 
superblock (ef = |ef B ). However one should note that this behaviour is not unusual in RG 
treatment of lattice systems. For example in standard RG at second order perturbation of 
free-fermion model 3 and isotropic Heisenberg model 14 , one obtains non-variational results. 
We have also obtained the ground state energy in the A = (xx model) case for other 
JlM - ratios, namely for the ratios t, tt and ^, the results are collected in table-2. It is 

n$B J 7' 9' 11 13' 

interesting to note that our results in these cases are above the exact values of ground state 
energy, which indicates that the non-variational result in the | case is due to the very small 
block and superblock sizes. Therfore we expect that for % and usb not very small, our 
method gives indeed variational results. Our data in table-2 show much accuracy compared 
with the standard RG results for any value of block sizes. These results show that the 
computed ground state energy indeed converges to the exact result faster than the standard 
RG method by increasing the ■^ B - ratios. 

° "SB 

3 Correlation Functions 

In this section we want to gain more insight on the effective Hamiltonian constructed by this 
modified RG. The xxz Hamiltonian in the massless region (0 < A < 1) has a non-degenerate 
ground state with a zero staggered magnetization in the thermodynamic limit N — >oo. The 



9 



staggered magnetization is defined as 

^H0|^E(-1)V*|0> (28) 

1=1 

In the RG formalizm the ground state |0) is replaced by T\0'), where |0') is the ground state 
of the effective Hamiltonian Hilbert space. Here we consider a chain of length 3 P and let p 
go to infinity. Then the staggered magnetization can be written in the following form 

m st = ^j:(-iy(0'\T^mO > ) (29) 

By dividing the chain into 3-sites blocks and consider T as the product of the embedding 
operator of each block, we can use eq (17) and write it again as 

where 

7i = 73 = (^) ; 72 = ( A2 )■ (30) 

Here we can continue this relation and replace |0') = T\0"), and so on. Then the staggered 
magnetization will be 

1 N 

^ t = ^E(-ir(0 (p) |7f^ (p) |0 (p) ) (31) 

which gives m st = 0, when p goes to infinity since (0 < 7« < 1) for (0 < A < 1). 

Another quantity which verifies the absence of long-range order in the massless region of 
the xxz model is the z-component of spin-spin correlation functions. It is defined as 

g(r) = <0|(-1)VM + ,|0> (32) 

for the anti-ferromagnetic ordering. In the classical Neel order, g(r) is a constant for any 
distance r, but quantum fluctuations can destroy Neel ordering. Despite the 2D Heisenberg 
model in the ID xxz model quantum fluctuations destroy the long-range order. 

We show that the constructed effective Hamiltonian will produce the above results. Again 
we assume a 3 P sites chain and let p — >oo. In this procedure sf is considered at the center 
of chain and the other blocks are expanded in both sides of it symmetrically. To calculate 
g(r), the ground state |0) is replaced by T|0') and eq.(30) is used for the effective operator 
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in the effective Hilbert space. The obtained g(r) is plotted in fig. 3a for xx model(A = 0) 
and compared with the exact results 15 . We have also plotted g(r) for A = 0.5 and A = 1 
in fig.3b. All of these graphs show good agreement with the exact results 15 . The good 
agreement of the correlation functions verify that the constructed effective Hamiltonian is a 
good approximate one for the low energy spectrum of the xxz model. 

4 Conclusions 

We have shown how by a simple modification of the standard RG method one can take care 
of the effect of boundary conditions in treating quantum lattice systems. We have applied 
this method to the spin \ xxz model in the massless region. We have shown our method 
give more accurate results than the standard RG for different value of block sizes. Our data 
of the ground state energy converge to the exact one faster than the standard RG results 
by increasing the block sizes. We have also calculated the staggered magnetizaton and the 
z-component of spin-spin correlation function in the massless region. These data confirm 
that we have obtained a good approximate Hamiltonian for the low energy spectrum which 
leads to no long-range order where, < A < I. These results are in good agreement with 
the exact results 11,12,15 . Our method although not as exact as the more sophisticated DMRG 
method has the merit that it can be implemented very easily on personal computers. One 
may hope that for a first estimate of many properties of these systems, our method or it's 
improvements (or adaptation for other lattice systems) may give acceptable results. This 
method can also be applied to other spin systems, to fermion models 10 , and more interestingly 
to 2-dimensional models (work on 2D anti- ferromagnetic Heisenberg model is in progress). 
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Tables 



Table-l. Ground state energy per site for standard RG, exact results and modified RG 
for xxz spin=| chain. 



A 


< E \ 

\ jjy ) standardRG 


( E \ 

\ jj\f J exact 


( E \ 

\ jjy jmodifiedRG 





-0.28284295 


-0.31830988 


-0.32529130 


0.1 


-0.29208390 


-0.32869220 


-0.33583683 


0.2 


-0.30163311 


-0.33956430 


-0.34603469 


0.3 


-0.31150403 


-0.35091277 


-0.35700253 


0.4 


-0.32171198 


-0.36272723 


-0.36813752 


0.5 


-0.33227452 


-0.37500000 


-0.37947617 


0.6 


-0.34321210 


-0.38772591 


-0.39120956 


0.7 


-0.35454900 


-0.40090224 


-0.40325429 


0.8 


-0.36631503 


-0.41452883 


-0.41546686 


0.9 


-0.37854863 


-0.42860840 


-0.42793815 


1 


-0.39130452 


-0.44314718 


-0.44076268 



Table-2. Ground state energy per site of standard RG, exact results and modified RG 
for A = at different value of block sizes. 



Exact result 


Standard RG 


Modified RG 


Type of SB and B 


-0.3183(^) 


n B = 3,-0.2828 


-0.3253 


usb = 5, n B = 3 




n B = 5, -0.2927 


-0.3113 


usb = 7, n B = 5 




n B = 7, -0.2983 


-0.3096 


usb = 9,n B = 7 




n B = 9,-0.3019 


-0.3098 


n S B = 11, n B = 9 




n B = 11,-0.3044 


-0.3103 


n S B = 13, n B = 11 
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Figure Captions 



Fig.l) Decomposition of the lattice to isolated blocks and considering neglected bonds 
as the effective interactions in the new Hilbert space. 
Fig. 2) Block and Super Block for the lattice chain. 

Fig. 3a) z-component of spin-spin correlation function g(r) = ((— l) r sfsf +r ) versus r for 
xx model(A = 0), modified RG and exact results. 

Fig. 3b) z-component of spin-spin correlation function g(r) = ((— l) r sf sf +r ) versus r for 
anisotropy parameter (A) equal to 0.5 and 1, obtained by modified RG. 
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